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Abstract The Effective One Body (EOB) formalism is an analytical approach which 
aims at providing an accurate description of the motion and radiation of coalescing 
binary black holes with arbitrary mass ratio. We review the basic elements of this 
formalism and discuss its aptitude at providing accurate template waveforms to be 
used for gravitational wave data analysis purposes. 



1 Introduction 

A network of ground-based interferometric gravitational wave (GW) detectors 
(LIGO/VIRGO/GEO/. . .) is currently taking data near its planned sensitivity [1]. 
Coalescing black hole binaries are among the most promising, and most exciting, 
GW sources for these detectors. In order to successfully detect GWs from coalesc- 
ing black hole binaries, and to be able to reliably measure the physical parameters 
of the source (masses, spins, . . .), it is necessary to know in advance the shape of 
the GW signals emitted by inspiralling and merging black holes. Indeed, the detec- 
tion and subsequent data analysis of GW signals is made by using a large bank of 
templates that accurately represent the GW waveforms emitted by the source. 

Here, we shall introduce the reader to one promising strategy toward having an 
accurate analytical' description of the motion and radiation of binary black holes, 
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which covers all its stages (inspiral, plunge, merger and ring-down): the Effective 
One Body approach [2, 3, 5, 4]. As early as 2000 [3] this method made several quan- 
titative and qualitative predictions concerning the dynamics of the coalescence, and 
the corresponding GW radiation, notably: (i) a blurred transition from inspiral to 
a 'plunge' that is just a smooth continuation of the inspiral, (ii) a sharp transition, 
around the merger of the black holes, between a continued inspiral and a ring-down 
signal, and (iii) estimates of the radiated energy and of the spin of the final black 
hole. In addition, the effects of the individual spins of the black holes were inves- 
tigated within the EOB [4, 6] and were shown to lead to a larger energy release 
for spins parallel to the orbital angular momentum, and to a dimensionless rotation 
parameter J/E^ always smaller than unity at the end of the inspiral (so that a Ken- 
black hole can form right after the inspiral phase). All those predictions have been 
broadly confirmed by the results of the recent numerical simulations performed by 
several independent groups [7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 
23, 24, 25, 26, 27, 28, 29] (for a review of numerical relativity results see also [30]). 
Note that, in spite of the high computer power used in these simulations, the cal- 
culation of one sufficiently long waveform (corresponding to specific values of the 
many continuous parameters describing the two arbitrary masses, the initial spin 
vectors, and other initial data) takes on the order of two weeks. This is a very strong 
argument for developing analytical models of waveforms. 

Those recent breakthroughs in numerical relativity (NR) open the possibility of 
comparing in detail the EOB description to NR results. This EOB/NR comparison 
has been initiated in several works [31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41]. The 
level of analytical/numerical agreement is unprecedented, compared to what has 
been previously achieved when comparing other types of analytical waveforms to 
numerical ones. In particular, Refs. [40, 41] have compared two different kind of 
analytical waveforms, computed within the EOB framework, to the most accurate 
gravitational waveform currently available from the Caltech-Cornell group, finding 
that the phase and amplitude differences are of the order of the numerical error. 

If the reader wishes to put the EOB results in contrast with other (Post-Newtonian 
or hybrid) approaches he can consult, e.g., [27, 28, 42, 43, 44, 45, 46, 47]. 

Before reviewing some of the technical aspects of the EOB method, let us indi- 
cate some of the historical roots of this method. First, we note that the EOB approach 
comprises three, rather separate, ingredients: 

1 . a description of the conservative (Hamiltonian) part of the dynamics of two black 
holes; 

2. an expression for the radiation-reaction part of the dynamics; 

3. a description of the GW waveform emitted by a coalescing binary system. 

For each one of these ingredients, the essential inputs that are used in EOB works 
are high-order post-Newtonian (PN) expanded results which have been obtained by 

to solve them. The important point is that, contrary to 3D numerical relativity simulations, nu- 
merically solving ODE's is extremely fast, and can therefore be done (possibly even in real 
time) for a dense sample of theoretical parameters, such as orbital (v = mimi/M,...) or spin 
(oi = Si/ Gm\ , 9i , (pi , . . .) parameters. 
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many years of work, by many researchers (see references below). However, one of 
the key ideas in the EOB philosophy is to avoid using PN results in their original 

"Taylor-expanded" form (i.e. cq + cj v + C2 + C3 H h c„ v"), but to use them 

instead in some resummed form {i.e. some non-polynomial function of v, defined 
so as to incorporate some of the expected non-perturbative features of the exact re- 
sult). The basic ideas and techniques for resumming each ingredient of the EOB are 
different and have different historical roots. Concerning the first ingredient, i.e. the 
EOB Hamiltonian, it was inspired by an approach to electromagnetically interacting 
quantum two-body systems introduced by Brezin, Itzykson and Zinn- Justin [48]. 

The resummation of the second ingredient, i.e. the EOB radiation-reaction force 
J^, was originally inspired by the Pade resummation of the flux function introduced 
by Damour, Iyer and Sathyaprakash [49]. Recently, a new and more sophisticated 
resummation technique for the radiation reaction force ^ has been introduced by 
Damour, Iyer and Nagar [50] and further employed in EOB/NR comparisons [40]. 
It will be discussed in detail below. 

As for the third ingredient, i.e. the EOB description of the waveform emitted by 
a coalescing black hole binary, it was mainly inspired by the work of Davis, Ruffini 
and Tiomno [51] which discovered the transition between the plunge signal and a 
ringing tail when a particle falls into a black hole. Additional motivation for the 
EOB treatment of the transition from plunge to ring-down came from work on the, 
so-called, "close limit approximation" [52]. 

Let us finally note that the EOB approach has been recently improved [37, 50, 40] 
by following a methodology consisting of studying, element by element, the physics 
behind each feature of the waveform, and on systematically comparing various 
EOB-based waveforms with 'exact' waveforms obtained by NR approaches. Among 
these 'exact' NR waveforms, it has been useful to consider the small-mass-ratio 
limit ^ V = mi m?/ (mi +m2)^ ^ 1, in which one can use the well controllable 'lab- 
oratory' of numerical simulations of test particles (with an added radiation-reaction 
force) moving in black hole backgrounds [35, 36]. 



2 IMotion and radiation of binary black holes: post- Newtonian 
expanded results 

Before discussing the various resummation techniques used in the EOB approach, 
let us briefly recall the 'Taylor-expanded' results that have been obtained by pushing 
to high accuracies the post-Newtonian (PN) methods. 

Concerning the orbital dynamics of compact binaries, we recall that the 2.5PN- 
accurate^ equations of motion have been derived in the 1980's [53, 54, 55, 56]. 



Beware that the fonts used in this chapter make the greek letter v (indicating the symmetric mass 
ratio) look very similar to the latin letter v indicating the velocity. 

' As usual ')!-PN accuracy' means that a result has been derived up to (and including) terms which 
ai'c ~ (v/c)-" ~ [GM/c^r)" fractionally smaller than the leading contribution. 
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Pushing the accuracy of the equations of motion to the 3PN (v/c)^) level proved 
to be a non-trivial task. At first, the representation of black holes by delta-function 
sources and the use of the (non diffeomorphism invariant) Hadamard regularization 
method led to ambiguities in the computation of the badly divergent integrals that 
enter the 3PN equations of motion [57, 58]. This problem was solved by using the 
(diffeomorphism invariant) dimensional regularization method {i.e. analytic contin- 
uation in the dimension of space d) which allowed one to complete the determina- 
tion of the 3PN-level equations of motion [59, 60]. They have also been derived by 
an Einstein-Infeld-Hoffmann-type surface-integral approach [61]. The 3.5PN terms 
in the equations of motion are also known [62, 63, 64]. 

Concerning the emission of gravitational radiation, two different gravitational- 
wave generation formalisms have been developed up to a high PN accuracy: (i) 
the Blanchet-Damour-Iyer formalism [65, 66, 67, 68, 69, 70, 71] combines a 
multipolar post-Minkowskian (MPM) expansion in the exterior zone with a post- 
Newtonian expansion in the near zone; while (ii) the Will-Wiseman-Pati formalism 
[72, 73, 74, 62] uses a direct integration of the relaxed Einstein equations. These for- 
malisms were used to compute increasingly accurate estimates of the gravitational 
waveforms emitted by inspiralling binaries. These estimates include both normal, 
near-zone generated post-Newtonian effects (at the IPN [66], 2PN [75, 76, 72], 
and 3PN [77, 78] levels), and more subtle, wave-zone generated (linear and non- 
linear) 'tail effects' [69, 79, 80, 71]. However, technical problems arose at the 3PN 
level. Similarly to what happened with the equation of motion, the representation 
of black holes by 'delta-function' sources causes the appearance of dangerously 
divergent integrals in the 3PN multipole moments. The use of Hadamard (par- 
tie finie) regularization did not allow one to unambiguously compute the needed 
3PN-accurate quadrupole moment. Only the use of the (formally) diffeomorphism- 
invariant dimensional regularization method allowed one to complete the 3PN-level 
gravitational-radiation formalism [82]. 

The works mentioned in this Section (see [83] for a detailed account and more 
references) finally lead to PN-expanded results for the motion and radiation of bi- 
nary black holes. For instance, the 3.5PN equations of motion are given in the form 
(a =1,2;/ =1,2,3) 



dt^ 



= A'=™'^+Af^, (1) 



where 



-Ao + c % + c (2) 



denotes the 'conservative' 3PN-accurate terms, while 

A'^'^ = c-^As+c-''At, (3) 

denotes the time-asymmetric contibutions, linked to 'radiation reaction' . 

On the other hand, if we consider for simplicity the inspiralling motion of a 
quasi-circular binary system, the essential quantity describing the emitted gravita- 
tional waveform is the phase of the quadrupolar gravitational wave amplitude 
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h{t) ~ fl(f)cos(0(f) + 5). PN theory allows one to derive several different func- 
tional expressions for the gravitational wave phase (j), as a function either of time 
or of the instantaneous frequency. For instance, as a function of time, admits the 
following explicit expansion in powers of = vc^ {tc~t)/ 5GM (where tc denotes a 
formal 'time of coalescence', M = m\ +m2 and v = mi m^jM^) 




with some numerical coefficients fl„,aj, which depend only on the dimensionless 
(symmetric) mass ratio v = m\m2/M^. The derivation of the 3.5PN-accurate ex- 
pansion (4) uses both the 3PN-accurate conservative acceleration (2) and a 3.5PN 
extension of the (fractionally) IPN-accurate radiation reaction acceleration (3) ob- 
tained by assuming a balance between the energy of the binary system and the 
gravitational-wave energy flux at infinity (see, e.g., [83]). 

Among the many other possible ways [84] of using PN-expanded results to pre- 
dict the GW phase ^{t), let us mention the semi-analytic T4 approximant [42, 32]. 
The GW phase defined by the T4 approximant happens to agree well during the 
inspiral with the NR phase in the equal mass case [27]. However, this agreement 
seems to be coincidental because the T4 phase exhibits significant disagreement 
with NR results for other mass ratios [39] (as well as for spinning black-holes [47]). 



3 Conservative dynamics of binary black holes: the Effective 
One Body approach 

The PN-expanded results briefly reviewed in the previous Section are expected to 
yield accurate descriptions of the motion and radiation of binary black holes only 
during their early inspiralling stage, i.e. as long as the PN expansion parameter 
Ye = GM/c^R (where R is the distance between the two black holes) stays signifi- 
cantly smaller than the value ^ ^ where the orbital motion is expected to become 
dynamically unstable ('last stable circular orbit' and beginning of a 'plunge' leading 
to the merger of the two black holes). One needs a better description of the motion 
and radiation to describe the late inspiral (say 7(, > -j^), as well as the subsequent 
plunge and merger. One possible strategy for having a complete description of the 
motion and radiation of binary black holes, covering all the stages (inspiral, plunge, 
merger, ring-down), would then be to try to 'stitch together' PN-expanded analytical 
results describing the early inspiral phase with 3D numerical results describing the 
end of the inspiral, the plunge, the merger and the ring-down of the final black hole, 
see, e.g.,Refs. [86, 32]. 

However, we wish to argue that the BOB approach makes a better use of all the 
analytical information contained in the PN-expanded results (l)-(3). The basic claim 
(first made in [2, 3]) is that the use of suitable resummation methods should allow 
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one to describe, by analytical tools, a sufficiently accurate approximation of the en- 
tire waveform, from inspiral to ring-down, including the non-perturbative plunge 
and merger phases. To reach such a goal, one needs to make use of several tools: (i) 
resummation methods, (ii) exploitation of the flexibility of analytical approaches, 
(iii) extraction of the non-perturbative information contained in various numerical 
simulations, (iv) qualitative understanding of the basic physical features which de- 
termine the waveform. 

Let us start by discussing the first tool used in the BOB approach: the systematic 
use of resummation methods. Essentially two resummation methods have been em- 
ployed (and combined) and some evidence has been given that they do significantly 
improve the convergence properties of PN expansions. The first method is the sys- 
tematic use of Pade approximants. It has been shown in Ref. [49] that near-diagonal 
Fade approximants of the radiation reaction force^ J? seemed to provide a good 
representation of ,^ down to the last stable orbit (which is expected to occur when 
R ^ 6GM/c^, i.e. when ~ g). In addition, a new route to the resummation of ^ 
has been proposed very recently in Ref. [50]. This approach, that will be discussed 
in detail below, is based on a new multiplicative decomposition of the metric mul- 
tipolar waveform (which is originally given as a standard PN series). In this case, 
Pade approximants prove to be useful to further improve the convergence properties 
of one particular factor of this multiplicative decomposition. 

The second resummation method is a novel approach to the dynamics of compact 
binaries, which constitutes the core of the Effective One Body (EOB) method. 

For simplicity of exposition, let us first explain the EOB method at the 2PN 
level. The starting point of the method is the 2PN-accurate Hamiltonian describing 
(in Arnowitt-Deser-Misner-type coordinates) the conservative, or time symmetric, 
part of the equations of motion (1) {i.e. the truncation A™"^ = Aq + c^^A2 + c^'^Aa, 
of Eq. (2)) say //2pn(^i ~ ^i^Pi^Pi)- going to the center of mass of the sys- 
tem (pi + P2= 0), one obtains a PN-expanded Hamiltonian describing the relative 
motion, q = - q^, p = = -pj. 

Hlft%9,P)=Ho{q,P) + ^H2iq,P) + ^H,iq,p), (5) 

where Ho{q,p) ~ + (with M = m\+ m2 and /i = mi OT2 /M) corresponds 
to the Newtonian approximation to the relative motion, while H2 describes IPN 
corrections and 2PN ones. It is well known that, at the Newtonian approximation, 
Ho{q,p) can be thought of as describing a 'test particle' of mass jj. orbiting around 
an 'external mass' GM. The EOB approach is a general relativistic generalization 
of this fact. It consists in looking for an 'external spacetime geometry' g^^y {x^ ; GM) 
such that the geodesic dynamics of a 'test particle' of mass within g^^l (x^ , GM) is 



* We henceforth denote by the Hamiltonian version of the radiation reaction term A'*'*, Eq. (3), 
in the (PN-expanded) equations of motion. It can be heuristically computed up to (absolute) 5.5PN 
[77, 81, 82] and even 6PN [85] order by assuming that the energy radiated in gravitational waves 
at infinity is balanced by a loss of the dynamical energy of the binary system. 



The Effective One Body description of the Two-Body problem 



7 



equivalent (when expanded in powers of 1 /c^) to the original, relative PN-expanded 
dynamics (5). 

Let us explain the idea, proposed in [2], for establishing a 'dictionary' between 
the real relative-motion dynamics, (5), and the dynamics of an 'effective' particle of 
mass jX moving in g^^l,{x^ ,GM). The idea consists in 'thinking quantum mechan- 
ically'^. Instead of thinking in terms of a classical Hamiltonian, H{q,p) (such as 
T/jpn'™' (5))' ^^'^ °f i'^s classical bound orbits, we can think in terms of the 
quantized energy levels E{n,£) of the quantum bound states of the Hamiltonian 
operator H{q,p). These energy levels will depend on two (integer valued) quan- 
tum numbers n and d. Here (for a spherically symmetric interaction, as appropriate 
to //■■'^lative-j^ £ parametrizes the total orbital angular momentum (L^ = £{£+ l)fP'), 
while n represents the 'principal quantum number' n = i + nr+l, where «r (the 'ra- 
dial quantum number') denotes the number of nodes in the radial wave function. The 
third 'magnetic quantum number' m (with —£< m < i) does not enter the energy 
levels because of the spherical symmetry of the two-body interaction (in the center 
of of mass frame). For instance, a non-relativistic Coulomb (or Newton!) interaction 

1 7 GMu 

gives rise to the well-known result 



E,M^-\,{^^)\ (7) 



which depends only on n (this is the famous Coulomb degeneracy). When consider- 
ing the PN corrections to Hq, as in Eq. (5), one gets a more compUcated expression 
of the form 



1 a 



2'"n2 



,2 



1 -L _ ("£11 I £22^ -X- _ ('£11 -u SZL -u £11 -u £12^ 



(8) 

where we have set a = GMfi /h = Gminii/Ti, and where we consider, for simphc- 
ity, the (quasi-classical) limit where n and £ are large numbers. The 2PN-accurate 
result (8) had been derived by Damour and Schafer [87] as early as 1988. The di- 
mensionless coefficients Cpq are functions of the symmetric mass ratio v = /z/M, 
for instance C40 = j (145 — 15v + V^). In classical mechanics {i.e. for large n and £), 
it is called the 'Delaunay Hamiltonian', i.e. the Hamiltonian expressed in terms of 
the action variables^ J = &i= 4^§ pcp d(p, and N = nh= I, + J, with /,- = § prdr. 

The energy levels (8) encode, in a gauge-invariant way, the 2PN-accurate relative 
dynamics of a 'real' binary. Let us now consider an auxiliary problem: the 'effec- 
tive' dynamics of one body, of mass pi, following a geodesic in some 'external' 



^ This is related to an idea emphasized many times by John Archibald Wheeler: quantum mechan- 
ics can often help us in going to the essence of classical mechanics. 

i„oc;^„ii„ a — E 
1 ■ 



^ We consider, for simplicity, 'equatorial' motions with m = i.e., classically, 9 
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(spherically symmetric) metric^ 

g^lJ^^dx^dx" = -A{R)c^ dT^ + B{R)dR^ + R^{de^ + sin^ e d(p^) . (9) 

Here, the a priori unknown metric functions A(R) and B{R) will be constructed in 
the form of expansions in GM/c^R: 

2 /^,^\ 3 



, , GM f GMV f GM 
, , , , GM , f GmV^ 



(10) 



where the dimensionless coefficients a„,b„ depend on v. From the Newtonian limit, 
it is clear that we should set ai = —2. By solving (by separation of variables) the 
'effective' Hamilton-Jacobi equation 

5eff = -4fff+/eff(P + 5eff(/?), (H) 

one can straightforwardly compute (in the quasi-classical, large quantum numbers 
limit) the Delaunay Hamiltonian f»'eff(A^eff,-4ff)j with A'eff = «eff^> -^eff = ^eff^ (where 
A^eff = JeS + lf, with if = 2i f pf dR, Pf = dS^f{{R) /dR). This yields a result of 
the form 



1 (X 

^eff(neff,4ff) ^ llC^ --^ ^ 

2 ntrr 



a2 / 



^eff 


^effx 






neff4ff 


«eff/ 


^eff 


^eff 


"^13 


'-22 


"eff^eff 


"eff*eff 



™4 / „eff „eff „eff „eff \ 

IX / Cio Lao Lti 



-31 _|_ ^ \ 

«eff4ff «eff/ 

(12) 

where the dimensionless coefficients c^^ are now functions of the unknown coeffi- 
cients an,bn entering the looked for 'external' metric coefficients (10). 

At this stage, one needs (as in the famous AdS/CFT correspondence) to define a 
'dictionary' between the real (relative) two-body dynamics, summarized in Eq. (8), 
and the effective one-body one, summarized in Eq. (12). As, on both sides, quantum 
mechanics tells us that the action variables are quantized in integers (A'reai = nh, 
Nes = n^ffh, etc.) it is most natural to identify « = «eff and £ = feff- One then still 
needs a rule for relating the two different energies E'^f^^"'^ and tfeff . Ref. [2] proposed 
to look for a general map between the real energy levels and the effective ones 
(which, as seen when comparing (8) and (12), cannot be directly identified because 



' It is convenient to write the 'external metric' in Schwarzschild-like coordinates. Note that the 
external radial coordinate R differs from the two-body ADM-coordinate relative distance R^'^^ — 
\q\. The transformation between the two coordinate systems has been determined in Refs. [2, 5]. 
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Mc 



real 
2 



lie 



■n+1, 



■n + lJ + l 



-n+ll 



n + 1, 



■ nl 



-nl 



Fig. 1 Sketch of the correspondence between the quantized energy levels of the real and effective 
conservative dynamics, n denotes the 'principal quantum number' (n = 11^ + 1.+ I, with ;v = 0, 1 , . . . 
denoting the number of nodes in the radial function), while (. denotes the (relative) orbital angular 



momentum (L 



l)/i ). Though the BOB method is purely classical, it is conceptually 



useful to think in terms of the underlying (Bohr-Sommeifeld) quantization conditions of the action 
variables and J to motivate the identification between n and ti in the two dynamics. 



they do not include the same rest-mass contribution ), namely 

/P / ^relative \ ^relative / prelative / prelative \ ^ \ 

(13) 

The 'correspondence' between the real and effective energy levels is illustrated in 
Fig. 1 

Finally, identifying S'ef[{n,£)/ jxc^ to f{E^f^'"^/ixc^) yields six equations, relat- 
ing the six coefficients c^^{a2,ay,bi,b2) to the six Cpg{v) and to the two energy 
coefficients ai and a2- It is natural to set bi = +2 (so that the linearized effective 
metric coincides with the linearized Schwarzschild metric with mass M = mi+ m2). 
One then finds that there exists a unique solution for the remaining five unknown 
coefficients 02,03,^2, Oli and a2- This solution is very simple: 

V 

fl2=0, a3=2v, ^2=4-6v, = a2=0- (14) 



^ Indeed E^°f = Mc^ + E'^^^"" = Mc- + Newtonian tenns + IPN/f^ + • • • , while (^effective = + 
A'+lPN/t- + ---. 
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Note, in particular, that the map between the two energies is simply 

^tc c-relative / crelative \ 2 4_ 2 4 

4ff ^^^g^/^^V^\ . m,c m c ^^^^ 



jlc^ jlc^ \ 2 pic^ J 2m\m2c' 

where s = = i^c^ + Elf^^^'^f is Mandelstam's invariant = -{px + ^2)'- 

Note also that, at 2PN accuracy, the crucial '^qq'' metric coefficient A (which 
fully encodes the energetics of circular orbits) is given by the remarkably simple 
PN expansion 

A2pn(/?) = 1 -2m + 2vm\ (16) 

where u = GM/{c'^R) and v = ll/M = m\ m2/{m\ +^2)^. 

The dimensionless parameter v = pi/M varies between (in the test mass 
limit mi <C 1112) and ^ (in the equal-mass case mi = m2). When v ^ 0, Eq. (16) 
yields back, as expected, the well-known Schwarzschild time-time metric coeffi- 
cient -goo'''* = 1 - 2i( = 1 - IGMjc^R. One therefore sees in Eq. (16) the role of v 
as a deformation parameter connecting a well-known test-mass result to a non trivial 
and new2PN result. It is also to be noted that the IPN EOB result Aipn(/?) = 1 - 2m 
happens to be v-independent, and therefore identical to A^*^*^* = 1 — 2m. This is re- 
markable in view of the many non-trivial v-dependent terms in the IPN relative 
dynamics. The physically real IPN v-dependence happens to be fully encoded in 
the function /(£) mapping the two energy spectra given in Eq. (15) above. 

Let us emphasize the remarkable simplicity of the 2PN result (16). The 2PN 
Hamiltonian (5) contains eleven rather compUcated v-dependent terms. After trans- 
formation to the EOB format, the dynamical information contained in these eleven 
coefficients gets condensed into the very simple additional contribution + 2vm^ in 
A(R\ together with an equally simple contribution in the radial metric coefficient: 
(A(/?)B(/?))2PN = 1 — 6 VM^. This condensation process is even more drastic when 
one goes to the next (conservative) post-Newtonian order; the 3PN level, i.e. ad- 
ditional terms of order ^(1/c^) in the Hamiltonian (5). As mentioned above, the 
complete obtention of the 3PN dynamics has represented quite a theoretical chal- 
lenge and the final, resulting Hamiltonian is quite complicated. Even after going to 
the center of mass frame, the 3PN additional contribution ■^H(,[q^p) to Eq. (5) in- 
troduces eleven new complicated v-dependent coefficients. After transformation to 
the EOB format [5], these eleven new coefficients get "condensed" into only three 
additional terms: (i) an additional contribution io A{R), (ii) an additional contribu- 
tion to B{R), and (iii) a ^{p'^) modification of the 'external' geodesic Hamiltonian. 
For instance, the crucial 3PN ^qq' metric coefficient becomes 

A3pn(/?) = 1 -2M + 2vM^ + fl4VM'*, (17) 

where M = GM/(c2;;), 

94 41 , 

04 = y-^7r^=^ 18.6879027, (18) 



The Effective One Body description of the Two-Body problem 



11 




Fig. 2 Various approximations and Fade resummation of the EOB radial potential A («), where u = 
GM/{c^R), for the equal-mass case V = 1/4. The vertical dashed lines indicate the corresponding 
(adiabatic) LSO location [2] defined by the condition d^S^ydR^ = dS^ydR = 0, where is the 
effective energy along the sequence of circular orbits (i.e., when — 0). 



while the additional contribution to B{R) gives 

D3pn(^) = (A(/?)B(/?))3pn = 1-6vm2 + 2(3v-26)vm^ . (19) 

Remarkably, it is found that the very simple 2PN energy map Eq. (15) does not need 
to be modified at the 3PN level. 

The fact that the 3PN coefficient 04 in the crucial 'effective radial potential' 
A3pn(7?), Eq. (17), is rather large and positive indicates that the v-dependent non- 
linear gravitational effects lead, for comparable masses (v ~ j), to a last stable (cir- 
cular) orbit (LSO) which has a higher frequency and a larger binding energy than 
what a naive scaling from the test-particle limit (v 0) would suggest. Actually, 
the PN-expanded form (17) of A3PN(i?) does not seem to be a good representation 
of the (unknown) exact function Aeob {R) when the (Schwarzschild-like) relative 
coordinate R becomes smaller than about 6GM/c^ (which is the radius of the LSO 
in the test-mass limit). In fact, by continuity with the test-mass case, one a priori 
expects that A3pn(/?) always exhibits a simple zero defining an EOB "effective hori- 
zon" that is smoothly connected to the Schwarzschild event horizon atR~ IGM/c^ 
when V ^ 0. However, the large value of the 04 coefficient does actually prevent 
A3PN to have this property when v is too large, and in particular when v = 1 /4, as it 
is visually explained in Fig. 2. The black curves in the figure represent the A func- 
tion at IPN (solid line), 2PN (dashed line) and 3PN (dash-dot fine) approximation: 
while the 2PN curve still has a simple zero, the 3PN does not, due to the large value 
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of 04. It was therefore suggested [5] to further resum^ A3pn(/?) by replacing it by a 
suitable Pade (P) approximant. For instance, the replacement of A3pn(/?) by 

A\iR) EE Pj [A3pn(/?)] = , 3 (20) 

1 +aiu + d2U +d^u^ 

ensures that the V = | case is smoothly connected with the v = limit, as Fig. 2 
clearly shows'". 

The use of Eq. (20) was suggested before one had any (reliable) non-perturbative 
information on the binding of close black hole binaries. Later, a comparison with 
some "waveless" numerical simulations of circular black hole binaries [89] has 
given some evidence that Eq. (20) is physically adequate. In Refs. [4, 89] it was also 
emphasized that, in principle, the comparison between numerical data and EOB- 
based predictions should allow one to determine the effect of the unknown higher 
PN contributions to Eq. (17). For instance, one can add a 4PN-like term +a5Vu^ 
or a 5PN-like term +a(,vu^ in Eq. (17), and then Pade the resulting radial function. 
The new resummed A potential will exhibit an explicit dependence on (at 4PN) 
or (flsjflg) (at 5PN), that is 



A\{R;ai,v)=Pl 



i3PN(^) + Vfl5M 



or 



A5(/?;fl5,fl6, V) A3pN(^) + Vfl5tr"' + Vfl6M 



5 , ,,„ ..6 



(21) 



(22) 



Comparing the predictions of A4 (/?; 05 , v) or Aj (7?; 05 , ag, v) to numerical data might 
then determine what is the physically preferred "effective" value of the unknown co- 
efficient 05 (if working at 4PN effective accuracy) or of the doublet {as,a^) (when 
including also 5PN corrections). For illustrative purposes. Fig. 2 shows the effect 
of the Pade resummation with = ag = and v = 1 /4. Note that the Pade re- 
summation procedure is injecting some "information" beyond that contained in the 
numerical values of the PN expansion coefficients fl„'s of A{R). As a consequence, 
the operation of Padeing and of restricting as and ag to the (3PN-compatible) val- 
ues as = = flg do not commute: A{(/?; 0,1/4) ^ A^(/?;0,0, 1/4) ^A\{R,\/A). 
In this respect, let us also mention that the 4PN 05 -dependent Pade approximant 
A\{R;as^v) exactly reduces to the 3PN Pade approximant A3 (/?; v) when 05 is re- 
placed by the following function of v 

3PN. . ^ v(3392- 123;r^)^ 
"5 18432(v-4) ■ ^^^^ 



^ The PN-expanded BOB building blocks A{R),B{R), . . . already represent a resummation of the 
PN dynamics in the sense that they have "condensed" the many terms of the original PN-expanded 
Hamiltonian within a very concise format. But one should not refrain to further resum the BOB 
building blocks themselves, if this is physically motivated. 

We recall that the coefficient ni and {di,d2,dj,) of the Pade approximant are determined by 
the condition that the first four tenns of the Taylor expansion of A3 in powers of w = GM/{c^R) 
coincide with A3PN . 
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Note that the value of the A3 -reproducing effective 4PN coefficient fl5™(v) in the 
equal mass case is ^^^(l /4) ~ —17.158031. This is numerically compatible with 
the value 05 = —17.16 quoted in Ref. [28] (but note that the correct A3 -reproducing 
4PN coefficient depends on the symmetric mass ratio v). Similarly, when work- 
ing at the 5PN level, A^^{R;as,a(,,v) exactly reduces to the 4PN Pade approximant 
A\{R;a5, v) when ae is replaced by the following function of both v and as: 

V (2304fl2 + 96 (3392 - 1237r2) as + (3776 - 123;r2) (32(3v + 94) - 1237r2)) 

24 [(3776 - 1237r2)v- 1536] ' 

(24) 

The use of numerical relativity data to constrain the values of the higher PN param- 
eters (05, ag) is an example of the useful flexibility [88] of analytical approaches; 
the fact that one can tap numerically-based, non-perturbative information to im- 
prove the EOB approach. The flexibility of the EOB approach related to the use of 
the as-dependent radial potential A\{R;as,v) has been exploited in several recent 
works [33, 37, 38, 39, 28, 41] focusing on the comparison of EOB-based wave- 
forms with waveforms computed via numerical relativity simulations. Collectively, 
all these studies have shown that it is possible to constrain as (together with other 
flexibility parameters related to the resummation of radiation reaction, see below) 
so as to yield an excellent agreement (at the level of the published numerical errors) 
between EOB and numerical relativity waveforms. The result, however, cannot be 
summarized by stating that as is constrained to be in the vicinity of a special nu- 
merical value. Rather, one finds a strong correlation between as and other parame- 
ters, notably the radiation reaction parameter Vpoie introduced below. More recently, 
Ref. [40] could get rid of the flexibility parameters (such as Vpoie) related to the 
resummation of radiation reaction, and has shown that one can get an excellent 
agreement with numerical relativity data by using only the flexibility in the doublet 
(asjag) (the other parameters being essentially fixed internally to the formalism). 
We shall discuss this result further in Sec. 5 below. 

The same kind of v-continuity argument discussed so far for the A function needs 
to be applied also to the D(/?)3pn function defined in Eq. (19). A straightforward 
way to ensure that the D function stays positive when R decreases (since it is D = 1 
when V 0) is to replace D3pn(^) by D^{R) = [D3pn(/?)], where indicates 
the (0,3) Pade approximant and explicitly reads 

^°(^)=l+6v«2_2(3v-26)v«3- (25) 

The resummation of A (via Pade approximants) is necessary for ensuring the exis- 
tence and v-continuity of a last stable orbit (see vertical lines in Fig. 2), as well as 
the existence and v-continuity of a last unstable orbit, i.e. of a v-deformed analog 
of the light ring R = 3GM/c^ when v ^ 0. We recall that, when v = 0, the Ught 
ring corresponds to the circular orbit of a massless particle, or of an extremely rel- 
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ativistic massive particle, and is technically defined by looking for the maximum 
of A{R)/R^, i.e. by solving {d/dR){A{R)/R^) = 0. When v 7^ and when consid- 
ering the quasi-circular plunge following the crossing of the last stable orbit, the 
"effective" meaning of the "v-deformed light ring" (technically defined by solving 
[d / dR){A{R : v)/R^) = 0) is to entail, in its vicinity, the existence of a maximum 
of the orbital frequency Q = d(p/dt (the resummation of D{R) plays a useful role 
in ensuring the v-continuity of this plunge behavior). 



4 Description of radiation-reaction effects in the Effective One 
Body approach 

In the previous Section we have described how the EOB method encodes the con- 
servative part of the relative orbital dynamics into the dynamics of an 'effective' 
particle. Let us now briefly discuss how to complete the EOB dynamics by defining 
some resummed expressions describing radiation reaction effects. One is interested 
in circularized binaries, which have lost their initial eccentricity under the influence 
of radiation reaction. For such systems, it is enough (as shown in [3]) to include a 
radiation reaction force in the equation of motion only. More precisely, we are 
using phase space variables r,pr,(p,p(p associated to polar coordinates (in the equa- 
torial plane 9 — Actually it is convenient to replace the radial momentum p,- by 
the momentum conjugate to the 'tortoise' radial coordinate /?* = / dR{B / A)^ 1^ , i.e. 
Pr^ = {A/Byi'^ Pr. The real EOB Hamiltonian is obtained by first solving Eq. (15) 
to get fi^g^f = in terms of (o'eff, and then by solving the effective Hamiltonian- 
Jacobi equation" to get (?eff in terms of the effective phase space coordinates q^^ 
and />gff. The result is given by two nested square roots (we henceforth set c = 1): 



with 73 = 2v (4 — 3 v). Here, we are using suitably rescaled dimensionless (effective) 
variables: r = R/GM, pr^, — Pr,/iJ., Pip — Pq,/ li GM, as well as a rescaled time t = 
T /GM. This leads to equations of motion (r, q),pr^,p(p) of the form 



" Completed by the ^(p*) terms that must be introduced at 3PN. 




(26) 



where 




(27) 
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d(p 
dt 

dr 
'dt 
dptp 
dt 

dpr, 

dt 



1/2 



<5Heob 



which expUcitly read 

Ap^ _ 



d(p 
dt 

dr 
dt 
dP(p 
dt 

dpr, 
dt 



Vr-HH^s 
B 



1 



vHH^ff 



Pr, 



2A , 



1/2 



1 



2v////eff 



-Z3 



A' 2A\ 4 



(28) 

(29) 
(30) 
(31) 

(32) 

(33) 
(34) 
(35) 



where A' = dA/dr. As explained above the EOB metric function A{r) is defined 
by Pade resumming the Taylor-expanded result (10) obtained from the matching 
between the real and effective energy levels (as we were mentioning, one uses a 
similar Pade resumming for D{r) = A{r)B{r)). One similarly needs to resum ^^f, 
i.e., the (j? component of the radiation reaction which has been introduced on the 
r.h.s. of Eq. (30). During the quasi-circular inspiral ^(f, is known (from the PN 
work mentioned in Section 2 above) in the form of a Taylor expansion of the form 



Taylor 



32 



(36) 



where V(f, = Q r^, and r^ = r[i/A(r,/7,|,)]'/^ is a modified EOB radius, with \jf being 
defined as 



. , 2 fdA{r) 



l+2v 



AW 1 + ^ -1 



(37) 



which generalizes the 2PN-accurate Eq. (22) of Ref. [90]. In Eq. (36) we have de- 
fined 



^Taylor^^,) = 1 + v) + v) A4(v) A5 (v) 

+A6(v,logy)v'^+A7(v)v^+A8(v = 0,logv) v*^ 



(38) 
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Fig. 3 The extreme-mass-ratio limit (v = 0): the Newton-normalized energy flux emitted by a 
particle on circular orbits. The figure illustrates the scattering of the standard Taylor expansion 
of the flux around the "exact" numerical result (computed up to = 6) obtained via perturbation 
theory. 



where we have added to the known 3.5PN-accurate comparable-mass result the 
small-mass-ratio 4PN contribution [91]. We recall that the small-mass contribution 
to the Newton-normalized flux is actually known up to 5.5PN order, i.e. to v'' in- 
cluded. The standard Taylor expansion of the flux, (38), has rather poor convergence 
properties when considered up to the LSO. This is illustrated in Fig. 3 in the small- 
mass limit V = 0. The convergence of the PN-expanded flux can be studied in detail 
in the v = limit, because in this case one can compute an "exact" result numeri- 
cally (using black hole perturbation theory [92, 93]). The "exact" energy flux shown 
in Fig. 3 is obtained as a sum over multipoles 

^-'^-^'^ £ L^im, (39) 
t=2 m=l 

where Ff,„ = Ff |,„| already denotes the sum of two equal contributions corresponding 
to -^m and — m (m 7^ as FjiQ vanishes for circular orbits). To be precise, the "exact" 
result exhibited in Fig. 3 is given by the rather accurate approximation F^^' obtained 
by choosing £niax = 6; i.e., by truncating the sum over t in Eq. (39) beyond ^ = 6. In 
addition, one normalizes the result onto the "Newtonian" (i.e., quadrupolar) result 
= 32/5(/x/M)^v"'. In other words, the solid line in Fig. 3 represents the quantity 

F = F^^yF^2- 

For clarity, we selected only three Taylor approximants; 3PN (v^), 3.5PN (v^) 
and 5.5PN (v'^). These three values suffice to illustrate the rather large scatter 
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among Taylor approximants, and the fact that, near the LSO, the convergence to- 
wards the exact value (solid line) is rather slow, and non monotonic. [See also Fig. 1 
in Ref. [94] and Fig. 3 of Ref. [49] for fuller illustrations of the scattered and non 
monotonic way in which successive Taylor expansions approach the numerical re- 
sult.] The results shown in Fig. 3 elucidate that the Taylor series (38) is inadequate 
to give a reliable representation of the energy loss during the plunge. That is the rea- 
son why the FOB formalism advocates the use of a "resummed" version of i.e. 
a nonpolynomial function replacing Eq. (38) at the rh.s. of the Hamilton's equation 
(and coinciding with it in in the vjc <C 1 limit). 

Two methods have been proposed to perform such a resummation. The first 
method, that strongly relies on the use of Pade approximants, was introduced 
by Damour, Iyer and Sathyaprakash [49] and, with different degrees of sophis- 
tication, has been widely used in the literature dealing with the FOB formal- 
ism [3, 6, 31, 35, 32, 33, 34, 36, 37, 38, 39, 28, 41]. The second resummation 
method has been recently introduced by Damour, Iyer and Nagar [50] and exploited 
to provide a self-consistent expression of the radiation reaction force in Ref. [40]. 
This latter resummation procedure is based on (i) a new multiplicative decomposi- 
tion of the gravitational metric waveform which yields a (ii) resummation of each 
multipolar contribution to the energy flux. The use of Pade approximants is a useful 
tool (but not the only one) that proves helpful to further improve the convergence 
properties of each multipolar contribution to the flux. The following two Sections 
are devoted to highlighting the main features of the two methods. For pedagogi- 
cal reasons the calculation is first done in the small-mass limit (v 0) and then 
generalized to the comparable mass case. 



4.1 Resummation ofF^^^^°'^ using a one-parameter family of Pade 
approximants: tuning Vpoie 



Following [49], one resums f' Taylor ^gjjjg j-jjg following Pade resummation ap- 
proach. First, one chooses a certain number Vpoie which is intended to represent the 
value of the orbital velocity at which the exact angular momentum flux would 
become infinite if one were to formally analytically continue J^^p along unstable 
circular orbits below the Last Stable Orbit (LSO): then, given Vpoie, one defines the 
resummed F{v^) as 



^resummed = ( 1 - 



-1 



p4 



^pole 



^T^=''°''(v^;v = 0) 



(40) 



where denotes a (4,4) Pade approximant. 

If one first follows the reasoning line of [49], and fixes the location of the pole in 

1 / a/3, one gets the 



the resummed flux at the standard Schwarzschild value v 



(v=0) _ 
pole 



result in Fig. 4. By comparison to Fig. 3, one can appreciate the significantly better 
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0.1 0.15 0.2 0.25 0.3 0.35 0.4 

V 



Fig. 4 The extreme-mass-ratio limit (v = 0). Pade resummation of the Taylor expandend energy 
flux of Fig. 3 as proposed in Ref. [49] with Vpoie = 1 / \/3. The sequence of Pade approximants is 
less scattered than the corresponding Taylor ones and closer to the exact result. 



(and monotonic) way in which successive Pade approximants approach (in La, norm 
on the full interval <x < jclso) the numerical result. Ref. [49] also showed that the 
observationally relevant overlaps (of both the "faithfulness" and the "effectualness" 
types) between analytical and numerical adiabatic signals were systematically better 
for Pade approximants than for Taylor ones. Note that this figure is slightly different 
from the corresponding results in panel (b) of Fig. 3 in [49] (in particular, the present 
result exhibits a better "convergence" of the v' ' curve). This difference is due to the 
new treatment of the logarithmic terms log a. Instead of factoring them out in 
front as proposed in [49], we consider them here (following [37]) as being part of 
the "Taylor coefficients" /„(logx) when Padeing the flux function. 

A remarkable improvement in the (Loo) closeness between /^Pade-resummed^^,^ 
pExact^^-^ can be obtained, as suggested by Damour and Nagar [37] (following ideas 
originally introduced in Ref. [97]), by suitably flexing the value of Vpoie- As pro- 
posed in Ref. [37], Vpoie is tuned until the difference between the resummed and the 
exact flux at the LSO is zero (or at least smaller than 10^^). The resulting closeness 
between the exact and tuned-resummed fluxes is illustrated in Fig. 5. It is so good 
(compared to the previous figures, where the differences were clearly visible) that 
we need to complement the figure with Table 1 . This table compares in a quantitative 
way the result of the "untuned" Pade resummation (vpoie = 1 / v^) of Ref. [49] to the 
result of the "vpoie-tuned" Pade resummation described here. Defining the function 
AF{v;vpoie) = Vpoie) — /'^"^'^'(v) measuring the difference between a 

resummed and the exact energy flux. Table 1 lists both the values of AF at v ~ v'lso 
and its Loo norm on the interval < v < vlso for both the untuned and tuned cases. 
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0.1 0.15 0.2 0.25 0.3 0.35 0.4 

V 



Fig. 5 The extreme mass ratio limit ( v = 0). Same of Fig. 4 hui flexing the value of the parameter 
Vpole so to improve the agreement with the exact result. 

Table 1 Errors in the flux of the two (untuned or tuned) Fade resummation procedures. From 
left to right, the columns report: the PN-order; the difference between the resummed and the exact 
flux, AF = - FExact^ 3j (jjg LSO, and the norm of AF, \ \AF\\^ (computed over the 

interval < v < I'lso). for I'pole ~ l/v'S; the flexed value of I'poie used here; AF at the LSO and 
the corresponding norm (over the same interval) for the flexed value of Vpoie . 

PN-order I IZiFllL^'^ Vp„.e Af;^^ \\AF\{i°'' 

3 (v") ^0048 0048 0.5334 7.06 x lO"'' 0.00426 

3.5 (v^) -0.051 0.051 0.5425 5.50x 10-"^ 0.00429 

5.5 (v") -0.022 0.022 0.5416 2.52x 10-^ 0.000854 



Note, in particular, how the Vpoie-flexing approach permits to reduce the Loo norm 
over this interval by more than an order of magnitude with respect to the untuned 
case. Note that the closeness between the tuned flux and the exact one is remarkably 
good (4.3 X IQ-^) already at the 3PN level. 

It has recently been shown in several works [37, 38, 39, 41] that Iht flexibility 
in the choice of Vpoie could be advantageously used to get a close agreement with 
NR data (at the level of the numerical error). We will not comment here any further 
on this paramete r-dependent resummation procedure of the energy flux and address 
the reader to the aforementioned references for further details. 
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4.2 Parameter-free resummation of waveform and energy jlux 

In this section we shall introduce the reader to the new resummation technique for 
the multipolar waveform (and thus for the energy flux) introduced in Ref. [36, 37] 
and perfected in [50]. The aim is to summarize here the main ideas discussed in [50] 
as well as to collect most of the relevant equations that are useful for implementa- 
tion in the EOB dynamics. To be precise, the new results discussed in Ref. [50] are 
twofold: on the one hand, that work generalized the £ = m = 2 resummed waveform 
of [36, 37] to higher multipoles by using the most accurate currently known PN- 
expanded results [99, 100, 101] as well as the higher PN terms which are known 
in the test-mass limit [95, 96]; on the other hand, it introduced a new resumma- 
tion procedure which consists in considering a new theoretical quantity, denoted as 
pim{x), which enters the {£,m) waveform (together with other building blocks, see 
below) only through its £-th power: hi,„ °^ (Pimix))^ ■ Here, and below, x denotes the 
invariant PN-ordering parameter x = (GM^2/c^)^/^. 

The main novelty introduced by Ref. [50] is to write the (£,m) multipolar wave- 
form emitted by a circular nonspinning compact binary as the product of several 
factors, namely 

f^l = ^n^^c,Uy)x'^^''''y'-'-- (f S^^T,y'-'pL (41) 

Here e denotes the parity of £ + m (e = n{i + m)), i.e. e = for "even-parity" (mass- 
generated) multipoles {£ + m even), and e = 1 for "odd-parity" (current-generated) 
ones {£ + m odd); n\J^ and c^+e ( v) are numerical coefficients; S^^ is a ^u-normahzed 
effective source (whose definition comes from the EOB formalism); T^,,, is a re- 
summed version [36, 37] of an infinite number of "leading logarithms" entering the 
tail ejfects [69, 103]; is a supplementary phase (which corrects the phase effects 
not included in the complex tail factor T^,,,), and, finally, {pimY denotes the £-th 
power of the quantity Pf m which is the new building block introduced in [50]. Note 
that in previous papers [36, 37] the quantity [pe,n) was denoted as fi,„ and we will 
mainly use this notation below. Before introducing explicitly the various elements 
entering the waveform (41) it is convenient to decompose he,„ as 

him = hf/^hf2, (42) 

where h^^^^ is the Newtonian contribution and /z^^^j = Sl,gT(,ne ft,„ represents a 

resummed version of all the PN corrections. The PN correcting factor n"^J^, as well 

as all its building blocks, has the structure h^^'>(,„ = 1 + ff[x). 

Entering now in the discussion of the explicit form of the elements entering 
Eq. (41), we have that the v-independent numerical coefficients are given by 
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.(0) ^ r....^^ 8^ / (£+l)(£ + 2) 
(2^+1)!! V £(£-1) 



- M %.. , ■ (43) 



(1) , 16;ri (2£+l)(£ + 2)(^^~m^) 

^ta (2^ + 1)!! y (2£-l)(£+ !)£(£- 1) ' ^ ^ 



while the v-dependent coefficients q_|_£(v) (such that |cf+£(v = 0)| = 1), can be 
expressed in terms of v (as in Ref. [99, 101]), although they are more conveniently 
written in terms of the two mass ratios X\ = m\/M and X2 = m2/M in the form 

= x|+'=-i + (-)'"xf+^-i. (45) 

In the second form of the equation we have used the fact that, as e = n{t + m), 
n{(. + e) = n{m). 

"(£) 

Let us turn now to discussing the structure of the S^g and T^,,, factors. To this aim, 
following Ref. [50], we recall that the along the sequence of BOB circular orbits, 
which are determined by the condition (9„ {A(m) [1 + Jq m^] } = 0, the effective BOB 
Hamiltonian (per unit jx mass) reads 



i/gg- = = ^A{u){\ + il m2) (circular orbits). (46) 

where the squared angular momentum is given by 

A'(u) 

il / 2AI ^^/ (circular orbits), (47) 

with the prime denoting d/du. Inserting this M-parametric representation of j' 
in Bq. (46) defines the M-parametric representation of the effective Hamiltonian 
H^ff{u). In the even-parity case (corresponding to mass moments), since the lead- 
ing order source of gravitational radiation is given by the energy density, Ref. [50] 
defined the even-parity "source factor" as 

s''^{x)=H,ff{x) e + m even, (48) 

where x = (GMQ/c^)'^^^. In the odd-parity case, they explored two, equally mo- 
tivated, possibilities. The first one consists simply in still factoring Hsff{x); i.e., in 
defining S^.^^ ~ Hes(x) also when i + m is odd. The second one consists in fac- 
toring the angular momentum ^ . Indeed, the angular momentum density EijkX^'fi^ 
enters as a factor in the (odd-parity) current moments, and ^ occurs (in the small-v 
limit) as a factor in the source of the Regge-Wheeler-Zerilli odd-parity multipoles. 
This leads us to define as second possibility 

S'^^g^'^ =]{x)=x'^l^j{x) i + m odd, (49) 
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where j denotes what can be called the "Newton-normalized" angular momentum, 
namely the ratio j{x) — j{x)/jN{x) with 7]v(x) = l/yx. In Ref. [50] the relative 
merits of the two possible choices were discussed. Although the analysis in the 
adiabatic v = limit showed that they are equivalent from the practical point of 
view (because they both yield waveforms that are very close to the exact numerical 
result) we prefer to consider only the 7-factorization in the following, that we will 
treat as our standard choice. 

The second building block in our factorized decomposition is the "tail factor" 
T(',„ (introduced in Refs. [36, 37]). As mentioned above, Ti,,, is a resummed version 
of an infinite number of "leading logarithms" entering the transfer function between 
the near-zone multipolar wave and the far-zone one, due to tail effects linked to 
its propagation in a Schwarzschild background of mass Madm = ^'eob- ^'■^ explicit 
expression reads 

_ r{£+l-2lk) I 2ihog(2kr„) ,.r.. 

hm- r(^+i) ' 

where ro = 2GM and k = GH^Q^mQ. and k = niQ. Note that k differs from A: by a 
rescaling involving the real (rather than the effective) EOB Hamiltonian, computed 
at this stage along the sequence of circular orbits. 

The tail factor Tg„ is a complex number which already takes into account some 
of the dephasing of the partial waves as they propagate out from the near zone to 
infinity. However, as the tail factor only takes into account the leading logarithms, 
one needs to correct it by a complementary dephasing term, e'^'^™, linked to sublead- 
ing logarithms and other effects. This subleading phase correction can be computed 
as being the phase 5e,„ of the complex ratio between the PN-expanded /z^J and the 
above defined source and tail factors. In the comparable-mass case (v ^ 0), the 3PN 
822 phase correction to the leading quadrupolar wave was originally computed in 
Ref. [37] (see also Ref. [36] for the v = limit). Full results for the subleading par- 
tial waves to the highest possible PN-accuracy by starting from the currently known 
3PN-accurate v-dependent waveform [101] have been obtained in [50]. 

The last factor in the multiplicative decomposition of the multipolar waveform 
can be computed as being the modulus f(,„ of the complex ratio between the PN- 
expanded /i^^^^ and the above defined source and tail factors. In the comparable mass 
case (v 7^ 0), the f22 modulus correction to the leading quadrupolar wave was com- 
puted in Ref. [37] (see also Ref. [36] for the v = limit). For the subleading partial 
waves, Ref. [50] explicitly computed the other /ft„'s to the highest possible PN- 
accuracy by starting from the currently known 3PN-accurate v-dependent wave- 
form [101]. In addition, as originally proposed in Ref. [37], to reach greater accu- 
racy the f(„,{x;vys extracted from the 3PN-accurate v 7^ results are completed 
by adding higher order contributions coming from the v = results [95, 96]. In 
the particular f22 case discussed in [37], this amounted to adding 4PN and 5PN 
V = terms. This "hybridization" procedure was then systematically pursued for 
all the other multipoles, using the 5.5PN accurate calculation of the multipolar de- 
composition of the gravitational wave energy flux of Refs. [95, 96]. Note that such 
hybridization procedure is not equivalent to the straightforward hybrid sum ansatz. 
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(v = 0) (where = h(,„/v) that one may have thought 



U„ — tknov/n(,,\ , ^higher 

to implement. 

In the even-parity case, the determination of the modulus ff„, is unique. In the 
odd-parity case, it depends on the choice of the source which, as explained above, 
can be connected either to the effective energy or to the angular momentum. We will 
consider both cases and distinguish them by adding either the label //or ^ to the 
corresponding ft,,,. Note, in passing, that, since in both cases the factorized effective 
source term (//eg- or ^) is a real quantity, the phases 5to,'s are the same. 

The above explained procedure defines the /fm's as Taylor-expanded PN series 
of the type 



fe„,{x-v) = 1 +c{<''"{v)x + 4'"{v)x^ 



„ftn 



{v ,\og{x))x^ + . . 



(51) 



Note that one of the virtues of our factorization is to have separated the half-integer 

is) 

powers of x appearing in the usual PN-expansion of /i^^J from the integer powers, 
the tail factor, together with the complementary phase factor e'^'"', having absorbed 
all the half-integer powers. In Ref. [39] all the /&„ 's (both for the H and choices) 
have been computed up to the highest available (v-dependent or not) PN accuracy. 
In the formulas for the /fo„'s given below we "hybridize" them by adding to the 
known v-dependent coefficients c{'"'(v) in Eq. (51) the v = value of the higher 
order coefficients: c^,"'{v = 0). The IPN-accurate f(„,'s for £ + m even and and also 
for i + m odd can be written down for all £. The complete result for the /ft„'s that 
are known with an accuracy higher than IPN are listed in Appendix B of Ref. [39]. 



Here, for illustrative purposes, we quote only the lowest and 
included. 



•odd, J 



up to £ = 3 



■ 114635 v^"* 227875 v2 



1 , , (2047v'- 6745 V -4288) 



41 



99792 
/ 36808 



9 

TT'V 



34625V 856 



V 2205 
/458816 
1^ 19845 



eulerlog2(jc) 



33264 96 

5391582359\ 4 
198648450 j ^ 



eulerlogofx) 

3696 105 ' 



eulerlog2 (x) 



93684531406 
893918025 



21428357\ 
727650 )' 



(52) 



J2d'-''>-^+[ 42 28j''+U52 126 9)'' 
/88404893 214 , A , 

+ (,11642400 -T05^"'^'"'°S'^"V" 

/6313 , , , , 33998136553\ 4 5, 
+ (^eulerlog,W-^^^^^jx^ + ^(x^), (53) 
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/33(x;v) = 1 



HH 



887v2 3401V 443 



330 



330 440 



53641811 



/ 147471561 78 , , , , , , , , , , ^ 
I 2802800 ' --l^ri°g3(-) j - + (^39 eulerlog3(x) - ^^^^^^ 



(54) 



, , 320v2-1115v + 328 39544v3 - 253768 v2+ 117215 v- 20496 , 



90(3v-l) 
110842222 104 , , , . > , 4. 



11880(3v-l) 



(55) 



/3i(a-;v) = 1 + 
/ 400427563 



2v 
T 
26 



13 
~6 



247 37 Iv 



198 



198 



V 75675600 21 



eulerlogj (x) I x + eulerlogj (x) — 



169 



1273 \ 2 
192 J"" 

12064573043\ 



63 



1816214400 /■ 



(56) 



For convenience and readability, we have introduced the following "eulerlog" func- 
tions eulerlog,„ (x) eulerlog,„ (x) = 7£ + log2 + -logx+logm, where Ye = 0.57721 . . . 
is Euler's constant. 

The decomposition of the total PN-correction factor /j^^^j into several factors is 
in itself a resummation procedure which has already improved the convergence of 
the PN series one has to deal with: indeed, one can see that the coefficients en- 
tering increasing powers of x in the /fo,'s tend to be systematically smaller than 
the coefficients appearing in the usual PN expansion of h^^^^. The reason for this is 
essentially twofold: (i) the factorization of T(„, has absorbed powers of mn which 
contributed to make large coefficients in hj,^J^, and (ii) the factorization of either Hes 
or j has (in the v = case) removed the presence of an inverse square-root sin- 
gularity located at x = 1/3 which caused the coefficient of x" in any PN-expanded 
quantity to grow as 3" as n To prevent some potential misunderstandings, 

let us emphasize that we are talking here about a singularity entering the analytic 
continuation (to larger values of x) of a mathematical function /i(x) defined (for 
small values of x) by considering the formal adiabatic circular limit. The point is 
that, in the v — *■ limit, the radius of convergence and therefore the growth with n 
of the PN coefficients of h{x) (Taylor-expanded at x = 0), are linked to the singu- 
larity of the analytically continued h{x) which is nearest to x = in the complex 
x-plane. In the v — *■ case, the nearest singularity in the complex x-plane comes 
from the source factor Hgf[{x) or j{x) in the waveform and is located at the light- 
ring xlr(v = 0) = 1/3. In the V ^0 case, the EOB formalism transforms the latter 
(inverse square-root) singularity in a more complicated ("branching") singularity 
where dHes/dx and dj/dx have inverse square-root singularities located at what 
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is called [3, 31, 33, 38, 37] the (Effective)'^ "EOB -light-ring", i.e., the (adiabatic) 

maximum of X2,xg^i;*(v) = [MQ^^^f^ > 1/3. 

Despite this improvement, the resulting "convergence" of the usual Taylor- 
expanded fem{xy& quoted above does not seem to be good enough, especially near 
or below the LSO, in view of the high-accuracy needed to define gravitational wave 
templates. For this reason, Refs. [36, 37] proposed to further resum the f22 (x) func- 
tion via a Pade (3,2) approximant, Pl{f22{x', v)}, so as to improve its behavior in the 
strong-field-fast-motion regime. Such a resummation gave an excellent agreement 
with numerically computed waveforms, near the end of the inspiral and during the 
beginning of the plunge, for different mass ratios [36, 38, 39]. As we were mention- 
ing above, a new route for resumming was explored in Ref. [50]. It is based on 
replacing /f,„ by its £-th root, say 

pim{x-v) = [fi„M-v)]'l^. (57) 

The basic motivation for replacing by p^,, is the following: the leading "Newtonian- 
level" contribution to the waveform /z^^J contains a factor co^r^^^v'^ where rharm is 
the harmonic radial coordinate used in the MPM formalism [66, 68] . When com- 
puting the FN expansion of this factor one has to insert the PN expansion of the 
(dimensionless) harmonic radial coordinate rhami, rhaim = x as 
a function of the gauge-independent frequency parameter x. The PN re-expansion of 
[niarm(j^)]^ then generates terms of the type x^' (1 +icix + ....). This is one (though 
not the only one) of the origins of IPN corrections in /zfo, and /&„ whose coefficients 
grow linearly with L The study of [50] has pointed out that these ^-growing terms 
are problematic for the accuracy of the PN-expansions. Our replacement of /ft,, by 
Piin is a cure for this problem. More explicitly, the the investigation of IPN correc- 
tions to GW amplitudes [66, 68, 99] has shown that, in the even-parity case (but see 
also Appendix A of Ref. [50] for the odd-parity case), 

1 ^ ' \ 3/ 2 2 q(v) q(v) q(v) 2(^+1)(2^ + 3)' 

(58) 

where C('(v) is defined in Eq. (45) and 

biiv)=xl+{-Yxl (59) 

Focusing on the v = case for simplicitly (since the v dependence of c{'"'(v) is 
quite mild [50]), the above result shows that the PN expansion of starts as 



Beware that this "Effective EOB-light-ring" occurs for a circular-orbit radius slightly larger 
than the purely dynamical (circular) EOB-light-ring (where H^fs and ^ would formally become 
infinite). 
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Fig. 6 Performance of the new resummation procedure described in Ref. [50]. The total GW flux 
F (up to £niax ~ 6) computed from inserting in Eq. (62) the factorized waveform (41) with the 
Taylor-expanded P(„,'s (with either 3PN or 5PN accuracy for P22) is compared with the "exact" 
numerical data. 



The crucial thing to note in this resuh is that as I gets large (keeping in mind that 
\m\ < i), the coefficient of x will be negative and will approximately range between 
—5£/4 and — £. This means that when i> 6 the IPN correction in f(,„ would by 
itself make fem{x) vanish before the (v = 0) LSO xlso =1/6. For example, for the 
£ = m = 6 mode, one has f^i^ix;Q) = 1 - 6.x(l + 1 1 /42) w 1 - 6x(l +0.26) which 
means a correction equal to —100% at .r = 1 /7.57 and larger than —100% at the 
LSO, namely f™il /6;0) « 1 - 1.26 = -0.26. This value is totally incompatible 
with the "exact" value /92^'^'(xlso) = 0.663 145 11 computed from numerical data in 
Ref. [50]. 

Finally, one uses the newly resummed multipolar waveforms (41) to define a 
resummation of the radiation reaction force is defined as 

^^^-^i^'^™-^', (61) 
where the (instantaneous, circular) GW flux Ft^max) jg defined as 

9 ^max 9 Anax ^ 

As an example of the performance of the new resummation procedure based on the 
decomposition of /if,,, given by Eq. (41), let us focus, as before, on the computation 
of the GW energy flux emitted by a test particle on circular orbits on Schwarzschild 
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spacetime. Figure 6 illustrates the remarkable improvement in the closeness between 
^New-resummed ^Exact rpjjg reader should Compare this result with the previous 
Fig. 3 (the straightforward Taylor approximants to the flux). Fig. 4 (the Pade re- 
summation with Vpoie = 1/V3) and Fig. 5 (the Vpoie-tuned Pade resummation). To 
be fully precise, Fig. 6 plots two examples of fluxes obtained from our new pim- 
representation for the individual multipolar waveforms /i^,,. These two examples 
differ in the choice of approximants for the l = m~2 partial wave. One example 
uses for P22 its 3PN Taylor expansion, Tt,[p22\, while the other one uses its 5PN 
Taylor expansion, Ts [p22\ ■ All the other partial waves are given by their maximum 
known Taylor expansion'^. Note that the fact that we use here for the Pfo,'s some 
straightforward Taylor expansions does not mean that this new procedure is not a 
resummation technique. Indeed, the defining resummation features of our procedure 
have four sources: (i) the factorization of the PN corrections to the waveforms into 
four different blocks, namely 5^^', r^,,,, e'^'"' and p^,^^ in Eq. (41); (ii) the fact the 

is by itself a resummed source whose PN expansion would contain an infinite 
number of terms; (iii) the fact that the tail factor is a closed form expression (see 
Eq. (50) above) whose PN expansion also contains an infinite number of terms and 
(iv) the fact that we have replaced the Taylor expansion of /f,„ = p^^^ by that of its 
£-th root, namely p^,. 

In conclusion, Eqs. (41) and (62) introduce a new recipe to resum the (v- 
dependent) GW energy flux that is alternative to the (vpoie-tuned) one given by 
Eq. (40). The two main advantages of the new resummation are: (i) it gives a better 
representation of the exact result in the v ^ limit (compare Fig. 6 to Fig. 5), and 
(ii) it is parameter-free: the only flexibility that one has in the definition of the wave- 
form and flux is the choice of the analytical representation of the function fii, like, 
for instance, P|{/22}, (Ts [p22])^, (Ts [p22])^, etc., (although Ref. [50] has pointed 
out the good consistency among all these choices). Note, that when V 7^ 0, the GW 
energy flux will depend on the choice of resummation of the radial potential A{R) 
through the Hamiltonian (for the even-parity modes) or the angular momentum (for 
the odd-parity modes). At the practical level, this means that the EOB model, im- 
plemented with the new resummation procedure of the energy flux (and waveform) 
described so far, will essentially only depend on the doublet of parameters (05,06), 
that can in principle be constrained by comparison with (accurate) numerical rela- 
tivity results. Contrary to the previous Vpoie -resummation of the radiation reaction, 
this route to resummation is free of radiation-reaction flexibility parameters. We 
will consider it as our "standard" route to the resummation of the energy flux in the 
following Sections discussing in details the properties of the EOB dynamics and 
waveforms. 



" We recall that Ref. [50] has also shown that the agreement improves even more when the Taylor 
expansion of the function P22 is further suitably Pade resummed. 
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5 Effective One Body dynamics and waveforms 

In this section we marry together all the EOB building blocks described in the pre- 
vious Sections and discuss the characteristic of the dynamics of the two black holes 
as provided by the EOB approach. In the following three subsections we discuss 
in some detail: (i) the set up of initial data for the EOB dynamics with negligible 
eccentricity (Sec. 5.1); (ii) the structure of the full Effective One Body waveform, 
covering inspiral, plunge, merger and ringdown, with the introduction of suitable 
Next-to-Quasi-Circular (NQC) effective corrections to it (and thus to the energy 
flux) (Sec. 5.2); (iii) the explicit structure of the EOB dynamics, discussing the so- 
lution of the dynamical equations. 



5.1 Post-post-circular initial data 

In this section we discuss in detail the so-called post-post-circular dynamical initial 
data (positions and momenta) as introduced in Sec. Ill B of [37]. This kind of (im- 
proved) construction is needed to have initial data with negligible eccentricity. Since 
the construction of the initial data is analytical, including the correction is useful to 
start the system relatively close and to avoid evolving the EOB equation of motion 
for a long time in order to make the system circularize itself. 

To explain the improved construction of initial data let us introduce a formal 
book-keeping parameter e (to be set to 1 at the end) in front of the radiation reaction 
in the EOB equations of motion. One can then show that the quasi-circular 
inspiralling solution of the EOB equations of motion formally satisfies 



Here, jo{r) is the usual circular approximation to the inspiralling angular momen- 
tum as explicitly given by Eq. (47) above. The order e ("post-circular") term Ui (r) is 
obtained by: (i) inserting the circular approximation p^ — jo{f) on the left-hand side 
(l.h.s) of Eq. (10) of [34], (ii) using the chain rule djo{r)/dt = {djo{r)/dr){dr/dt), 
(iii) replacing dr/dt by the right-hand side (rh.s) of Eq. (9) of [34] and (iv) solving 
for pr, at the first order in e. This leads to an explicit result of the form (using the 
notation defined in Ref. [34]) 



where the subscript indicates that the rh.s. is evaluated at the leading circular 
approximation e ^ 0. The post-circular EOB approximation {jo,ni) was intro- 
duced in Ref. [3] and then used in most of the subsequent EOB papers [6, 31, 



P^^Hr) + e^j2ir)+0{e'), 
Pr, = eni{r) + e^7i3{r) + 0{e^). 



(63) 



(64) 




(65) 
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32, 33, 34, 35]. The post-post-circular approximation (order £^), introduced in 
Ref. [37] and then used systematically in Ref. [38, 39, 40], consists of: (i) for- 
mally solving Eq. (35) with respect to the explicit appearing on the r.h.s., (ii) 
replacing p^* by its post-circular approximation, Eq. (65), (iii) using the chain rule 
dn\{r)/dt = {dni{r) / dr){dr / dt), and (iv) replacing dr/dt in terms of n\ (to lead- 
ing order) by using Eq. (33). The result yields an explicit expression of the type 
pi ~ jl{r)[\+£^k2{r)] of which one finally takes the square root. In principle, this 
procedure can be iterated to get initial data at any order in e. As it will be shown 
below, the post-post-circular initial data {jo E^ki, K\ ) are sufficient to lead to 
negligible eccentricity when starting the integration of the EOB equations of motion 
at radius r = R/{GM) = 15. 



5.2 Effective One Body waveforms 

At this stage we have essentially discussed all the elements that are needed to com- 
pute the EOB dynamics obtained by solving the EOB equation of motion, Eqs. (32)- 
(35). The dynamics of the system yields a trajectory {q{t),p{t)) = {(p{t),r{t),p^{t),pi-, 
in phase space. The (multipolar) metric waveform during the inspiral and plunge 
phase, up to the EOB "merger time" f,„ (that is defined as the maximum of the orbital 
frequency Q,) is a function of this trajectory, i.e. /j^^^^P'""^'^ = /j'"^'''""^'^ (?),/?(?)). 
Focusing only on the dominant £ = m ^2 waveform, the waveform that describes 
the full process of the binary black hole coalescence (i.e., inspiral, plunge, merger 
and ringdown) can be split in two parts: 

• The insplunge waveform: /i'"'*P'™se(f)^ computed along the EOB dynamics up to 
merger, which includes (i) the resummation of the "tail" terms described above 
and (ii) some effective parametrization of Next-to-Quasi-Circular effects. The 
€ = m = 2 metric waveform explicitly reads 

(^) ^T'^'^t) = v4°)c2(v)./';22(V;x)/2^Q^y2,-2 g^^) (66) 

where the argument x is taken to be (following [90]) x = = [raQ.)^ (where 
Tft) was introduced in Eq. (36) above). The resummed version of /22 enter- 
ing in h22{x) used here is given by the following Pade- re summed function 
72^2^ = ^2 [/^^'™(jc; v)]. In the waveform h22 above we have introduced (follow- 
ing [40]) a new ingredient, a "Next-to-Quasi-Circular" (NQC) correction factor 
of the form''* ^ 

/^2^<=(fll,fl2) = l+fll7^^+fl2-^, (67) 



Note that one could also similarly improve the subleading higher-multipolar-order contributions 
to J?^. In addition, other (similar) expressions of the NQC factors can be found in the literature [38, 
39,41]. 
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where and 02 free parameters that have to be fixed. A crucial facet of the 
new EOB formaHsm presented here consists in trying to be as predictive as pos- 
sible by reducing to an absolute minimum the number of "flexibility parameters" 
entering our theoretical framework. One can achieve this aim by "analytically" 
determining the two parameters 01,02 entering (via the NQC factor Eq. (67)) the 
(asymptotic) quadrupolar EOB waveform ^/if^^ (where R = R/M) by imposing: 
(a) that the modulus reaches, at the EOB-determined "merger time" ?,„, 

a local maximum, and (b) that the value of this maximum EOB modulus is equal 
to a certain (dimensionless) function of v, (p(v). In Ref. [40] we calibrated ^(v) 
(independently of the EOB formalism) by extracting from the best current Nu- 
merical Relativity simulations the maximum value of the modulus of the Numer- 
ical Relativity quadrupolar metric waveform I^/ijt'I. Using the data reported in 
[29] and [39], and considering the "Zerilli-normalized" asymptotic metric wave- 
form ^22 Rh22/V24, we found ^(v) ~ 0.3215v(l - 0.131(1 -4v)). Our re- 
quirements (a) and (b) impose, for any given A (m) potential, two constraints on 
the two parameters a 1,02- We can solve these two constraints (by an iteration 
procedure) and thereby uniquely determine the values of a 1,02 corresponding 
to any given A{u) potential. In particular, in the case considered here where 
A{u) = A{u;a5,a(,,v) this uniquely determines ai,a2 in function of as,a^ and 
V. Note that this is done while also consistently using the "improved" version of 
/122 given by Eq. (66) to compute the radiation reaction force via Eq. (62). 
• a simplified representation of the transition between plunge and ring-down 
by smoothly matching (following Refs. [36]), on a (2/7+ l)-toothed "comb" 
,t,n — 5,t,n,tm + 5, . . . ,fm + p5) Centered around a matching time f„,, 
the inspiral-plus-plunge waveform to a ring-down waveform, made of the super- 
position of several'^ quasi-normal-mode complex frequencies, 

(i^) = e-<^'-'-^ , (68) 

with (7^ = ttiv + icON, and where the label refers to indices {£,i\m,n), with 
{i,m) — (2,2) being the Schwarzschild-background multipolarity of the consid- 
ered (metric) waveform /2f,„, with n = 0,1,2... being the 'overtone number' 
of the considered Kerr-background Quasi-Normal-Mode, and £' the degree of 
its associated spheroidal harmonics S(i,„{aa ,0). As discussed in [3] and [36], 
and already mentioned above, the physics of the transition between plunge and 
ring-down (which was first understood in the classic work of Davis, Ruffini and 
Tiomno [51]) suggests to choose as matching time f„,, in the comparable-mass 
case, the EOB time when the EOB orbital frequency i2(f) reaches its maximum 
value. 

Finally, one defines a complete, quasi-analytical EOB waveform (covering the 
full process from inspiral to ring-down) as: 



" Refs. [36, 38] use p = 2, i.e. a 5-teethed comb, and, con'espondingly, 5 positive-frequency Kerr 
Quasi-Normal Modes. 
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hTit) = 0(?™-O/4"(O + 0(f-f«.)/^22' (0, (69) 

where 0(f) denotes Heaviside's step function. The final resuh is a waveform that 
only depends on the two parameters (as, as) which parametrize some flexibility on 
the Pade resummation of the basic radial potential A (m), connected to the yet uncal- 
culated (4PN, 5PN and) higher PN contributions. 



5.3 Effective One Body dynamics 

We conclude this section by discussing the features of the typical EOB dynamics ob- 
tained by solving the EOB equation of motion Eqs. (32)-(35) with post-post-circular 
initial data. The resummation of the radiation reaction force uses the multiplicative 
decomposition of hf,^ given by Eq. (41) with NQC correction to the £ ~ m ~ 2 
multipole given by Eq. (67). We fix the free parameters to the model to be — 0, 
ag = —20 (see below why) while ai and 02 are obtained consistently according to 
the iteration procedure discussed above. The system is started at ro = 15 and (po = 0. 
The post-post-circular initial data give /7° = 4.31509298 and = -0.00109847. 
The result of the outcome of the integration of the EOB equation of motion is dis- 
played in Fig. 7 together with the trajectory (top-left panel) and the orbital frequency 
(bottom-right panel). On this plot we remark two things. First, the fact that the or- 
bital frequency has a maximum at time f,„ = 3522 that identifies, in EOB, the merger 
(and matching) time. Second, the fact that p,-^ tends to a finite value after the merger 
(contrary to pr, that would diverge), yielding a more controllable numerical treat- 
ment of the late part of the EOB dynamics. 



6 Effective One Body and Numerical Relativity waveforms 

So far we have seen that (at least) two different EOB models (of dynamics and 
waveforms) are available. They differ, essentially, in the way the resummation of the 
GW energy flux yielding the radiation reaction force is performed. The first EOB 
model, that we will refer to as the "old" one, basically uses a Pade-resummation 
of the energy flux with an external parameter Vpoie that must be fixed in some way. 
The second EOB model, that we will refer to as the "improved" one, uses a more 
sophisticated resummation procedure of the energy flux, multipole by multipole, 
in such a way that the final result depends explicitly only on the same parameters 
{as , fle) that are used to parametrize higher PN contribution to the conservative part 
of the dynamics. 

In the last three years, the power of the "old" EOB model has been exploited in 
various comparisons with numerical relativity data, aiming at constraining in some 
way the space of the EOB flexibility parameters (notably represented by and 
Vpoie) by looking at regions in the parameter space where the agreement between the 
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Fig. 7 BOB dynamics for 05 = and = —20. Clockwise from the top left panel,the panels 
report: the trajectory, the radial separation r(f), the radial momentum Pr, (conjugate to r*), the 
orbital frequency £2{t), the angular momentum P(p{t) and the orbital phase (p{t). 



numerical and analytical waveforms is at the level of numerical error. For example, 
after a preliminary comparison done in Ref. [31], Buonanno et al. [33] compared 
restricted BOB waveforms"^ to NR waveforms computed by the NASA-Goddard 
group, showing that it is possible to tune the value of so as to have a good 
agreement between the two set of data. In particular, for 05 ~ 60 and Vpoie given 
according to the (nowadays outdated) suggestion of Ref. [49], in the equal-mass 
case (v = 1/4), they found that the dephasing between (restricted) EOB and NR 
waveforms (covering late inspiral, merger and ring-down) stayed within ±0.030 



The terminology "restricted" refers to a waveform which uses only the leading Newtonian ap- 
proximation, h''^f \ to the wavefonn 
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GW cycles over 14 GW cycles. In the case of a mass ratio 4 : 1 (v = 0.16), the 
dephasing stayed within ±0.035 GW cycles over 9 GW cycles. 

Later, the resummed factorized BOB waveform of Eq. (66) above within the "old" 
EOB model has been compared to several set of equal-mass and unequal-mass NR 
waveforms: (i) in the comparison with the very accurate inspiralling simulation of 
the Caltech-Cornell group [27] the dephasing stayed smaller than ±0.001 GW cy- 
cles over 30 GW cycles (and the amplitudes agreed at the ^ 10^^ level) [37]; (ii) 
in the comparison [38] with a late-inspiral-merger-ringdown NR waveform com- 
puted by the Albert Einstein Institute group, the dephasing stayed smaller than 
±0.005 GW cycles over 12 GW cycles; (iii) in the (joint) comparison [39] between 
EOB and very accurate equal-mass inspiralling simulation of the Caltech-Cornell 
group [27] and late-inspiral-merger-ringdown waveform for 1:1, 2:1 and 4:1 mass 
ratio data computed by the Jena group it was possible to tune the EOB flexibility 
parameters (notably a=, and Vpoie) so that the dephasing stayed at the level of the 
numerical error The same "old" model, with resummed factorized waveform, and 
the parameter-dependent (using Vpoie) resummation of radiation reaction force, was 
recently extended by adding 6 more flexibility parameters to the ones already intr- 
duced in Refs. [37, 39], and was "calibrated" on the high-accuracy Caltech-Cornell 
equal-mass data [41]. This calibration showed that only 5 flexibility parameters (as, 
Vpoie and three parameters related to non-quasi-circular corrections to the waveform 
amplitude) actually suffice to make the "old" EOB and NR waveform agree, both 
in amplitude and phase, at the level of the numerical error (this multi-flexed EOB 
model brings in an improvement with respect to the one of Refs. [37, 39] espe- 
cially for what concerns the agreement between the waveform amplitude around the 
merger). 

Recently, Ref. [40] has introduced and fully exploited the possibilities of the "im- 
proved" EOB formalism described above, taking advantage of: (i) the multiplica- 
tive decomposition of the (resummed) multipolar waveform advocated in Eq. (41) 
above, (ii) the effect of the NQC corrections to the waveform (and energy flux) 
given by Eq. (66), and, most importantly, (iii) the parameter-free resummation of 
radiation reaction In Ref. [40] the (asjflg) -dependent predictions made by 
the "improved" formalism were compared to the high-accuracy waveform from 
an equal-mass BBH (v = 1/4) computed by the Caltech-Cornell group [29], (and 
now made available on the web). It was found that there is a strong degeneracy 
between as and ag in the sense that there is an excellent EOB-NR agreement for 
an extended region in the (fl5,fl6)-plane. More precisely, the phase difference be- 
tween the EOB (metric) waveform and the Caltech-Cornell one, considered between 
GW frequencies Mcol = 0.047 andMojR = 0.31 (i.e., the last 16 GW cycles before 
merger), stays smaller than 0.02 radians within a long and thin banana-like region 
in the (ci5,ci6)-plane. This "good region" approximately extends between the points 
(flSjfle) = (0,-20) and {a$,a^) = (—36, +520). As an example (which actually 
lies on the boundary of the "good region"), we have followed [40] in considering 
here the specific values 05 = 0,fl6 = —20 (to which correspond, when v = 1/4, 
fli = -0.036347,(32 = 1.2468). We henceforth use M as time unit. 
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Fig. 8 This figure illustrates the comparison between the "improved" EOB wavefonn (quadrupolar 
(£ = m = 2) metric waveform (66) with parameter-free radiation reaction (61) and with 05 = O.og = 
—20) with the most accurate numerical relativity waveform (equal-mass case) nowadays available. 
The phase difference between the two is 40 < ±0.01 radians during the entire inspiral and plunge. 
Ref. [40] has shown that this agreement is at the level of the numerical error. 



This result relies on the proper comparison between NR and EOB time series, 
which is a delicate subject. In fact, to compare the NR and EOB phase time-series 
(l>ri' {(nr) and ^^2^{tEOB) one needs to shift, by additive constants, both one of the 
time variables, and one of the phases. In other words, we need to determine T and a 
such that the "shifted" EOB quantities 

fEOB=fEOB + T, 0;f°«- 02^^ + (70) 

"best fit" the NR ones. One convenient way to do so is first to "pinch" the EOB/NR 
phase difference at two different instants (corresponding to two different frequen- 
cies). More precisely, one can choose two NR times f|^'^,f^'^, which determine two 
corresponding GW frequencies'^ coi = (B^2^(ff"^), CO2 = tO^'^(f^'^), and then find 
the time shift z{(Oi , (O2) such that the shifted EOB phase difference, between (Oi and 

O}2,A0^''^{X) = (/);f OB (^EOB ) - (/EOB ) = ,^|OB EOB ^ ^) _ ^EOB EOB ^ ^) 

is equal to the corresponding (unshifted) NR phase difference Zi^"^*^ = 0|2^(f2"^) — 
(^22^(f j''^). This yields one equation for one unknown (t), and (uniquely) determines 
a value z{(Oi ,(02) of t. [Note that the ©2 ^ fi)i = C0,„ limit of this procedure yields 
the one-frequency matching procedure used in [27].] After having so determined T, 
one can uniquely define a corresponding best-fit phase shift a(a)i , fih) by requiring 

that, say, ^^^^(fjEOB) = ,^EOB(fjEOB) ^ ^ ^ <^)2NR(fNR), 



' Alternatively, one can start by giving oneself coi , CO2 and determine the NR instants fj , at 
which they are reached. 
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Having so related the EOB time and phase variables to the NR ones we can 
straigthforwardly compare all the EOB time series to their NR correspondants. In 
particular, we can compute the (shifted) EOB-NR phase difference 

Figure 8 compares'^ (the real part of) our analytical metric quadrupolar waveform 
"^22°^/^ to the corresponding (Caltech-Cornell) NR metric waveform ^^^/ V. This 
NR metric waveform has been obtained by a double time-integration (following the 
procedure of Ref. [39]) from the original, publicly available, curvature waveform 
Such a curvature waveform has been extrapolated both in resolution and in 
extraction radius. The agreement between the analytical prediction and the NR result 
is striking, even around the merger See Fig. 9 which closes up on the merger The 
vertical line indicates the location of the EOB-merger time, i.e., the location of the 
maximum of the orbital frequency. 

The phasing agreement between the waveforms is excellent over the full time 
span of the simulation (which covers 32 cycles of inspiral and about 6 cycles of 
ringdown), while the modulus agreement is excellent over the full span, apart from 
two cycles after merger where one can notice a difference. More precisely, the phase 



° The two frequencies used for this comparison, by means of the "two-frequency pinching tech- 
nique" mentioned above, are Mffli = 0.047 and M^y, = 0.31. 
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Fig. 10 Comparison between Numerical Relativity and EOB metric waveform for the 2: 1 mass 
ratio. 



difference, = — ^mSric remains remarkably small (~ ±0.02 radians) dur- 

ing the entire inspiral and plunge (0)2 = 0.31 being quite near the merger). By com- 
parison, the root-sum of the various numerical errors on the phase (numerical trun- 
cation, outer boundary, extrapolation to infinity) is about 0.023 radians during the 
inspiral [29]. At the merger, and during the ringdown, A(j) takes somewhat larger 
values (~ ±0. 1 radians), but it oscillates around zero, so that, on average, it stays 
very well in phase with the NR waveform (as is clear on Fig. 9). By comparison, 
we note that [29] mentions that the phase error linked to the extrapolation to infinity 
doubles during ringdown. We then note that the total "two-sigma" NR error level 
estimated in [29] rises to 0.05 radians during ringdown, which is comparable to the 
EOB-NR phase disagreement. In addition, Ref. [40] compared the "improved" EOB 
waveform to accurate numerical relativity data (obtained by the Jena group [39]) on 
the coalescence of unequal mass-ratio black-hole binaries. Fig. 10 shows the re- 
sult of the EOB/NR waveform comparison for a 2:1 mass ratio, corresponding to 
V = 2/9. When 05 = 0, ag = -20 one finds ai ^ -0.017017 and 02 = 1.1906. 
Again, the agreement is excellent, and within the numerical error bars. 

Finally, Ref. [40] explored another aspect of the physical soundness of the EOB 
analytical formalism: the triple comparison between (i) the NR GW energy flux at 
infinity (which was computed in [28]); (ii) the corresponding analytically predicted 
GW energy flux at infinity (computed by summing |/ift„P over ); and (iii) (mi- 
nus) the mechanical energy loss of the system, as predicted by the general EOB 
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Fig. 11 The triple comparison between Numerical Relativity and EOB GW energy fluxes and the 
EOB mechanical energy loss. 



formalism, i.e. the "work" done by the radiation reaction ^mechanical — 

Q^^. This 

comparison is shown in Fig. 11, which should be compared to Fig. 9 of [28]. We 
kept here the same vertical scale as [28] which compared the NR flux to older ver- 
sions of (resummed and non-resummed) analytical fluxes and needed such a ±10% 
vertical scale to accomodate all the models they considered. [The horizontal axis 
is the frequency (D of the differentiated metric waveform h22-] By contrast, we see 
again the striking closeness (at the ~ 2 x 10^^ level) between the EOB and NR GW 
fluxes. As both fluxes include higher multipoles than the (2,2) one, this closeness 
is a further test of the agreement between the improved EOB formalism and NR re- 
sults. [We think that the ~ 2a difference between the (coinciding) analytical curves 
and the NR one on the left of the Figure is due to uncertainties in the flux computa- 
tion of [28], possibly related to the method used there of computing h.] Note that the 
rather close agreement between the analytical energy flux and the mechanical en- 
ergy loss during late inspiral is not required by physics (because of the well-known 
"Schott term" [104]), but is rather an indication that h(,„ can be well approximated 
by —imQh{,„ 
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7 Conclusions 

We have reviewed the basic elements of the Effective One Body (EOB) formaHsm. 
This formaHsm is still under development. The various existing versions of the EOB 
formalism have all shown their capability to reproduce within numerical errors the 
currently most accurate numerical relativity simulations of coalescing binary black 
holes. These versions differ in the number of free theoretical parameters. Recently 
a new "improved" version of the formalism has been defined which contains essen- 
tially only two free theoretical parameters. 
Among the successes of the EOB formalism let us mention: 

1. An analytical understanding of the non-adiabatic late-inspiral dynamics and of 
its "blurred" transition to a quasi-circular plunge; 

2. The surprising possibility to analytically describe the merger of two black holes 
by a seemingly coarse approximation consisting of matching a continued inspiral 
to a ringdown signal; 

3. The capability, after using suitable resummation methods, to reproduce with 
exquisite accuracy both the phase and the amplitude of the gravitational wave 
signal emitted during the entire coalescence process, from early-inspiral, to late- 
inspiral, plunge, merger and ringdown; 

4. The gravitational wave energy flux predicted by the EOB formalism agrees, 
within numerical errors, with the most accurate numerical-relativity energy flux; 

5. The ability to correctly estimate (within a 2% error) the final spin and mass of 
nonspinning coalescing black hole binaries [this issue has not been discussed in 
this review, but see Ref. [34]]. 

We anticipate that the EOB formalism will also be able to provide an accurate 
description of more complicated systems than the nonspinning BBH discussed in 
this review. On the one hand, we think that the recently improved EOB frame- 
work can be extended to the description of (nearly circularized) spinning black 
hole systems by suitably incorporating both the PN-expanded knowledge of spin 
effects [105, 106, 108] and their possible EOB resummation [4, 107]. On the other 
hand, the EOB formalism can also be extended to the description of binary neutron 
stars or mixed binary systems made of a black hole and a neutron star [109, 110]. 
An important input for this extension is the use of the relativistic tidal properties of 
neutron stars [111, 112, 113] 

Finally, we think that the EOB formalism has opened the realistic possibility of 
constructing (with minimal computational resources) a very accurate, large bank of 
gravitational wave templates, thereby helping in both detecting and analyzing the 
signals emitted by inspiralling and coalescing binary black holes. Though we have 
had in mind in this review essentially ground-based detectors, we think that the 
EOB method can also be applied to space-based ones, i.e., to (possibly eccentric) 
large mass ratio systems. 
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